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Abstract 

We extend previous work on applying the e expansion to universal properties of a cold, dilute 
Fermi gas in the unitary regime of infinite scattering length. We compute the ratio ^ = /^/ep of 
chemical potential to ideal gas Fermi energy to next-to-next-to-leading order (NNLO) in e = 4 — d, 
where d is the number of spatial dimensions. We also explore the nature of corrections from the 
order after NNLO. 
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I. INTRODUCTION AND RESULTS 



For a number of years, it has been a challenge to compute the properties of a dilute 
Fermi gas with infinite scattering length [1]. This is known as the unitary regime and is 
relevant to cold systems for which the interparticle separation is very small compared to 
the two-particle scattering length a, but very large compared to any other distance scales 
r characterizing two-body interactions: a 3> n~^^^ ^ r, where n is the density. In the 
unitary regime, the only dimensionful parameter is the number density n, and so all physical 
quantities in this limit should be determined by dimensional analysis and universal constants 
of proportionality. Dimcnsionless ratios will be universal [2]. There has been much interest 
in recent years in attempting to compute such universal constants of the unitary regime (see, 
e.g., Refs. [3, 4] and references therein). 

The problem is nonperturbative in three spatial dimensions. However, inspired by earlier 
work of Nussinov and Nussinov [5] on the behavior of the unitary regime as a function 
of spatial dimension d, it was recently realized [6] that, with an appropriately formulated 
perturbation theory, a perturbative solution is possible in d=A—e spatial dimensions when 
e <^ 1. Results can be expressed as an asymptotic series in e, analogous to the e expansion 
methods that have been used with great success for 30 years to determine critical exponents 
in a variety of second-order phase transitions. One may then extrapolate to the case of three 
dimensions, e = 1. For the case of cold, dilute Fermi gases at infinite scattering length, it is 
found [6] that 

i = l^ = le'l''^ [£3/2 _ 0.0492 e^/^ + 0{e"'')\ 

^ 1^3/2 + Xe5/2 In g _ 0.0246 e'/^ + • • • , (1.1) 

- = 2 + 0(e), (1.2) 
A* 

A 2 

_ = - -0.691 + 0(e), (1.3) 

jJi C 

where is the chemical potential, A is the gap for fermionic excitations, and Eq is the value 
of /2m for fermionic excitations with the minimum energy A. The Fermi energy ep is 
defined as the Fermi energy of an ideal Fermi gas with the same density n as the strongly 
interacting gas under consideration. 

The ratio ^ can be equivalently expressed as an energy density ratio^ 

i-y. (1-4) 

where £ and £q are the energy densities in the interacting and noninteracting cases, respec- 
tively, at equal number density n. For experimental relevance, ^ can also be expressed as 
[7, 8] 



^rel \ 2 



Ef 



(harmonic trap) (1.5) 



^ This can be proven simply by using thermodynamics and scahng at unitarity [2]. 
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for a system in a harmonic trap (in the hmit of an arbitrarily wide trap), where the "release" 
energies i?''*^^ and Eq*^^ are the total system energies in the interacting and noninteracting cases 
at equal total particle number A^. 

Historically, in the simplest application of the e expansion, one computes the first two or 
three terms in the e expansion and sees if they are reasonably well behaved for e—1, which 
corresponds to three spatial dimensions. Using this method, the above expansions would 
give the estimates [6] 

£- 0.475, —^2, — -1.31. (1.6) 

In this method, it is important to quit when one is ahead: Because the e expansion is 
asymptotic, higher-order terms eventually grow. The simple procedure is to stop including 
higher-order terms when this happens. In more sophisticated applications of the e expansion, 
however, critical exponents have been determined fairly precisely for some phase transitions 
by combining high-order e expansions, information about the large-order asymptotic behav- 
ior, and knowledge of the behavior at or near lower dimensions such as d=2, to fit results as a 
function of dimension using Borel-Pade approximations [9-12].^ A similar procedure has re- 
cently been carried out for the ratio ^ — /i/ep in Ref. [13] using next-to- leading-order results 
for the e expansion about both four and two spatial dimensions (except that the high-order 
asymptotic behavior of the e expansion about four dimensions is currently unknown) . 

The goal of the current paper is to take the step of computing the next term in the e 
expansion about four spatial dimensions for £ at zero temperature, and to learn something 
about the analytic structure in e by showing that new, non-trivial logarithms of e appear at 
yet higher orders. We find 

£ = [e3/2 _ 0.04916 e^/^ - 0.95961 e'^^ - | e^/^ ]^, + o{e'/^)] . (1.7) 

One may expand 

e^/^d ^ I + lelne + e^^ln\ + j^lne) + 0{eHn^ e) (1.8) 

if desired, as in Eq. (1.1). If one intends to naively set e=l, this expansion is unnecessary 
since Ine then vanishes and e'^/^'' = 1 order by order in e. We will comment on the large 
relative size of the next-to-next-to-leading-order (NNLO) correction at the end of the paper 
and discuss there the implications of our result for Borel-Pade extrapolations to d—3. 

In the remainder of this paper, we explain our calculation of C, to NNLO in the e expansion. 
In the next section, we briefly review the formalism and diagrammatic rules developed in Ref. 
[6] for applying the e expansion to this problem. At the beginning of Sec. Ill, we display 
the diagrams that need to be evaluated to push the calculation of the effective potential 
(which is later used to determine to NNLO in e. Because the efficient evaluation of some 
diagrams is challenging, we will then take the time to explain our methods in detail. Results 
for all the diagrams are summarized in Appendix A. In Sec. IV, we put everything together 
to determine £ as in Eq. (1.7). Finally, in Sec. V, we explain how the diagrammatic e power 
counting of Ref. [6] would break down, due to infrared issues, if wc proceeded to yet one 
higher order in e than the calculation reported in this paper. We then show how the proper 
power counting of e can be restored. Finally, we discuss the implications of our result to 
extrapolating the value of £ to d—S in Sec. VI. 



^ For a textbook overview, see chapters 28 and 41 of Ref. [14]. 
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II. REVIEW 



We will generally follow the conventions and diagrammatic methods of Ref. [6], which we 
review here. One starts with the Lagrangian 



£ = + ^ ^ + fji^^as^ - —(f)*(f) + ¥a+^(f) + (2.1 

\ 2m J Co 



where is a Hubbard- Stratonovich field, ^ = , V'j)^ is a two-component Nambu-Gor'kov 
field, a± = |((Ti =t io"2), and cti 2,3 are the Pauli matrices. The constant Co determines the 
scattering length. In three spatial dimensions, 

V/^^. (2.2) 



Ana Co J (27r)^ p 

The integral on the right hand side of Eq. (2.2) is ultraviolet divergent (and is so in any 
dimension d > 2). In a physical system there is always an upper momentum cutoff, for 
example, the inverse of the range of potential. However, if the system is insensitive to the 
cutoff, then cq should always appear in physical observables in the same combination as in 
Eq. (2.2), so that observables can be expressed in terms of the scattering length. 

Technically, it is cumbersome to carry the integral in Eq. (2.2) across our formulas, so we 
will use dimensional regularization. The technical advantage of this regularization scheme 
is that the integral fd'^p/p'^ vanishes in any dimension, and so the connection between cq 
and a becomes very simple. In particular, infinite scattering length corresponds to Cq = 00, 
so that there is no term quadratic in in Eq. (2.1). As long as the physics is insensitive to 
short distances, dimensional regularization will give the same result as other regularization 
schemes such momentum cutoff. 

Following Ref. [6], we expand 4> about its superfluid expectation value {4>) = 4>o as 

= 0o + fl'</?, (2.3) 

where g is chosen to give the dynamics of 99 a conventional normalization at leading order 
in e. In particular, if one computes the small-momentum expansion of the (p self-energy H 
of Fig. 1 to leading order in e, one finds 

H(po,p) = H(0) + ^ (^-po + 1^) [1 + 0{e)] + 0{pI,p'). (2.4) 

The choice g'^m'^ = 87r^e[l + 0(e)] will then make the momentum dependence above into a 
conventionally normalized kinetic term for a nonrelativistic particle of mass = 2m for 
d — A. Reference [6] found it convenient to define 



(2.5) 



2 _ Stt e / m(j)Q 



m^ \ 2n 



^ Sec Ref. [15] for a standard textbook treatment of dimensional regularization, with emphasis in high- 
energy physics. For condensed matter applications see, for instance, the textbook treatment in Ref. [14]. 
For a few examples of use in the theory of cold, dilute atomic gases, see Ref. [16]. 
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FIG. 1: The one- loop scalar self-energy H, where solid lines represent Nambu-Gor'kov fermion 
fields * and dashed lines represent the scalar field variation ip. (There is a similar diagram where 
one of the external scalar arrows is reversed, mixing with <^*, which we will discuss later and 
denote n.) 

Reference [6] then reorganized the Lagrangian in the case of infinite scattering length as 
a sum £ = £o + -^1 + -^2 corresponding to an unperturbed Lagrangian of a free fermion 
field ^ and a free scalar field 93, plus perturbations L\ -|- L^- 

£0 = (idt + ^ + cr+00 + c^-0o) ^ + + ^) 9?, (2.6a) 

= 5f^V+*(^ /x^Vg* 2/x(^V, (2.6b) 

£2 = ( — — 1inp*i^. (2.6c) 

Here, L\ can be thought of as the interaction terms. As explained in Ref. [6], £2 should be 
employed as counterterms to the one-loop diagrams shown in Fig. 2.^ Feynman rules are 
shown in Fig. 3, where 

2 

no=-po + ^. (2.7) 

(This is our one deviation from the notation of Ref. [6]. We find it convenient to define our 
rio as the negative of their Ho. We have added the hat over 11 to distinguish it and avoid 
notational confusion on this point.) The unperturbed propagators generated by are 

^^'^ 00 PO + 6p - isj - vl-El^ie \ -00 PO - 6 J ^2-^^ 

and 

^(po,p) = 1^-— (2-9) 

Po - 2^P + 



^ The split-up (2.6) might appear unconventional, but the basic idea behind it. as in many oth(T occasions 
in condensed matter physics, is to resum divergent graphs. The graphs that needs to be resunimed in our 
case are multiple insertions of the fermion loop (Fig. 1) into the •p propagator. One could, in principle, 
formulate a set of Feynman rules where the y propagator is the inverse of the fermion loop, and fermion 
loop insertion into the scalar propagator is forbidden by hand. However, the resulting propagator would 
be a very complicated function of momentum and chemical potential. For practical calculations it is much 
more efficient to give the ip propagator a simpler form, equal to the inverse of the leading 1/e piece of the 
fermion loop [Eq. (2.4)], and have it corrected in higher loops. To do that at the formal level, the split-up 
(2.6) is introduced. 
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(a) --> -( )-->-- + ->-■->- 



o 



(b) .->..( ).->.. + ..>..0..>. 



FIG. 2: The rule for using combining divergent subdiagrams with counterterms from £2 to achieve 
a simple perturbative expansion in e [6] . 



vertices 



) ' ) ) ' ) ->*->- 

A 

igG+ igG- i Uq 

) X ) ->-x->-- 

j|j.a3 2i\i -2i\i 



propagators 



> 



iG iD 
FIG. 3: Feynman rules from Eqs. (2.6) [6]. 

where 

ep^l^, E^^{el + <t>lf'\ (2.10) 

By analyzing the effective potential l^(0o) for the expectation (pQ, the minimum is found at 
00 ~ Ai/Cj or equivalently ~ e0o- The correspondence between the diagrammatic expansion 
and the e expansion can then be codified by treating 0o as 0{e'^) and each insertion of fi or 
[see Eq. (2.5)] as 0(e). With one exception, this identification gives the relative importance 
of each diagram of the effective potential to NNLO in the e expansion (the order relevant 
for the current calculation), provided one uses counterterms according to Fig. 2. The one 
exception is the one- loop fermion diagram with a single /i insertion, shown in Fig. 4(b), which 
produces a 1/e that compensates for the 11 and which is not canceled by any counterterm 
diagram. 

A word of review is in order regarding the nature of the problem above and below four 
dimensions. In = 4 — e, the corrections to mean-field theory for the unitary Fermi gas are 
controlled by powers of e. For d > 4, on the other hand, the problem becomes qualitatively 
different, as the short distance scale R (the range of the potential) renders the problem non- 
universal, and ill defined as i? — > 0. This is easy to see if one solves the two-body problem in 
arbitrary dimension, as done by Nussinov and Nussinov [5]. [Precisely at c? = 4, mean-field 
theory acquires logarithmically small corrections of order 1/ log(I//i?), where L ~ n~^/^ is 
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(a) A (b) ^ (c) v> 

O OO 

FIG. 4: One-loop diagrams through O(e^). Scalar loops arc not shown smcc these vanish due to 
the retarded nature of the propagators We will use the notation (a) V^'^ , (b) V^^^ , and (c) y'^^^^ . 







FIG. 5: Two-loop diagrams. The counter-term diagrams for (a-c) — a single scalar loop with an 
appropriate £2 counterterm — are not shown because they vanish due to the retarded nature of the 
scalar propagator. Our notation is (a) v!f^ ^ and (b-d) v!f^ . 



the typical particle separation.^] 



^ One can understand this logarithm in the spirit of Nussinov and Nussinov by recalling that the unitary 
limit corresponds to the prosonco of a zero-energy bound state, whose wave function will be ^E* (x r^^^ 
outside of the range R of the potential. The normalization integral / d'^r|\I'p ~ / d'^r r'^^'^'^ for the 
total probability is UV convergent in = 4 — e dimensions but UV logarithmically divergent in = 4 
dimensions, where it introduces a logarithmic dependence on R. More technically, if one follows the 
d = 4 — e derivations of Nishida and Son, briefly reviewed here, the source of the small parameters ~ e 
and /i/00 ~ e of the expansion about mean-field theory come from logarithmically divergent (in d = 4) 
integrals: respectively, the 1/e in the self-energy (2.4) above [from Fig. 1] and in the potential (4.1) [from 
Fig. 4(b)]. These integrals / rf'^p/p^ are momentum-space versions of / d'^rr^"^''. Imagine roughly cutting 
ofl[ the integrals in the ultraviolet at the scale r ~ i? and p where the effective theory breaks down 

and one would need a treatment of the details of the two-body potential. In the calculations of Nishida 
and Son and in this paper, the infrared is cut off by the distance scale s ~ {rn<^{^)~^l'^ associated with the 
condensate (^q. In d = 4, we then see that the role of 1/e is replaced by ln(s/i?). (The resulting solution 
for 00 will relate s and L by a power of this logarithm, and so \r\.{sjR) ^ ln{L/R) up to corrections 
suppressed by inverse powers of the logarithm.) We also learn that, in d = 4 — e dimensions with e small, 
R must be exponentially tiny in order to be in the universal regime. Specifically, in order for R not to 
significantly affect / d'^p/p'^, one must have R <C e~^/'^s ~ e^/^e"^/*^!/. 
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-L -K' 




-L -K> 




-L -K' 




(h) 

r 



FIG. 6: Three-loop diagrams (including counterterms to bosonic self-energies), (h) vanishes due to 
the retarded nature of the scalar propagator, but we find it useful to include to maJce^more obvious 
the cancellation of UV divergences. Our notation is (a) V^^^^^ (b-d) V^^'^\ (e) V^^'^'^^ and (f-h) 



(nn) 



(a) 



.-<■-. 



(b) 



(c) 



*" -X- 



FIG. 7: Examples of diagrams which vanish because of the purely retarded nature of the scalar 
propagator D{P). 



III. EVALUATING THE NEXT-TO-NEXT-TO-LEADING-ORDER POTENTIAL 

A. The diagrams 

Figures 4-6 show the nontrivial diagrams that determine the effective potential V{(f)o) 
through NNLO in e, together with our conventions for labeling momenta. In the figure 
captions, we explain the notation we will use for the contribution of various classes of these 
diagrams to the effective potential. We have generally not included various one-loop scalar 
diagrams which vanish simply because of the retarded nature of the scalar propagator D of 
Eq. (2.9), such as those shown in Fig. 7.^ 



^ Naively, these diagrams can be seen to vanish by closing the loop frequency integration in the upper half 
plane, which contains no poles. There is a technical caveat, however, in that the contribution from the 
semicircle at infinity cannot be ignored in all cases. This can give (/)o-independent contributions to the 
effective potential which vanish in dimensional regularization and which in other schemes correspond to 
operator ordering issues, such as whether /x multiplies the Wigner-ordered number operator ^{a^a + aa^) 
(corresponding to a naive application of the Feynman rules) or the correct normal-ordered operator a^a. 
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B. Cross diagram 



As our first example, we will begin by discussing how to efficiently evaluate the contribu- 
tion Vi®^(0o) of the cross diagram of Fig. 6(a) to the full effective potential V{(po) at O(e^). 
This is a three-loop diagram and its leading contribution is already O(e^) because of the four 
explicit factors of the coupling g. We may therefore ignore any other e dependence at the 
order of interest and so can evaluate the loop integrals for exactly d=4. (For this diagram, 
the loop integrals converge for (i=4, which we will see explicitly.) Our basic approach will 
be to evaluate all of the frequency integrals. One can then scale all dimensionful parameters 
out of the remaining momentum integrals. We then perform the dimensionless momentum 
integrals numerically. 

Figure 6(a) gives 

-il/g^®) = -(ig)^ f iGn{P) iG22{Q) iG2i{-K) iGi2{P+Q+K) iD{P+K) iD{-Q-K), 

JpQK 

(3.1) 

where the overall minus sign on the right-hand side is for the fermion loop. We will use 
capital letters P to stand for (po, p), where po is frequency and p is spatial momentum, and 
we will use the short-hand notations 

f ^ fdpoj^ [ ^ fd^ 

Jp "-J 27r (27r)<^""' J^"-J {27ry"'- ^"''^ 
Since G21 — G12, this is the same as 

) ^ ^ /■ G'n(P) G22{Q) G^i-K) G,2{P+Q+K) D{P+K) D{-Q-K). (3.3) 

^ JpQK 



1. The frequency integrals 



One could simply use the expressions (2.8) and (2.9) for the propagators and now do the 
three frequency integrals (po) Qo-, ^0) by brute force. This approach is tedious and yields 
complicated expressions with many terms requiring significant effort to simplify. It also 
naturally produces terms with energy denominators such as {Eq — Ek + |eq+fe)~^ which look 
like they prodiicc singularities for certain momenta (e.g., q= — k in this example), but all 
such singularities turn out to cancel between different terms in the final result. 

There is, however, a method for carrying out the frequency integration which directly 
produces much tidier results. The first step is to rewrite 



Gu{P) = 

G22{Q) = 

Gu{S) = 



Po + ep 



1 



pl-E^ + ze 2Ep 

go - ^ 1 

qi-El + ie~ 2Eq 

-00 _ ^00 

sl-El + ie~ 2E, 



E„ + e. 



+ 



E - € 



{po - Ep + ie) (pi 



Ep — ie) 



Eq - eq 



+ 



E„ 



(qo-Eq + ie) {qo + Eq - ie) 
1 1 



(so - Es + ie) (so + Eg- ie) 



(3.4a) 
(3.4b) 
(3.4c) 



i.e., to decompose the propagators into retarded and advanced parts. The integrand in Eq. 
(3.3) then sphts into 2^ = 16 different terms. We shall see that many of these terms trivially 
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p 



+ - 
^ + <: 



FIG. 8: Split of fermion propagators into retarded and advanced terms according to Eq. (3.4). 




(a) (b) 



FIG. 9: (a) One of the 16 terms generated from Fig. 6a by the expansions (3.4). (b) A vanishing 
loop in this diagram. 

vanish, and others are related by symmetry. It is useful to give a graphical depiction of 
these different terms by schematically rewriting Eqs. (3.4) as in Fig. 8. The first and second 
terms on the right of this figure denotes the {po — Ep + ie)^^ and (po — Ep — ie)"^ terms on 
the right-hand side of Eqs. (3.4). The + and — signs in Fig. 8 denote the sign of ie, and the 
direction of the arrows on the right-hand side of Fig. 8 correspondingly represent the flow 
of time (forward for a retarded propagator, backward for an advanced one). 

In this new notation. Fig. 9(a) shows an example of one of the 2^^ terms contained in 
the original diagram of Fig. 6(a). It is easy to see that this term vanishes, because there 
exists a loop. Fig. 9(b), where all the arrows have the same orientation. If we do that loop 
integration first, then we can close it in a half plane where there are no poles, and we obtain 
zero. There are only six terms that do not contain a similar vanishing loop, and they are 
shown in Fig. 10. The terms represented by the bottom row are related to those of the top 
row by the change of variables 

{P,Q,K)^-{Q,P,K), (3.5) 

which is a symmetry of the original integrand (3.3).'' This just represents a change of 
integration variables, and so the contribution of the second row to the potential will equal 
that of the first row. We therefore need only evaluate three terms, corresponding to Figs. 
10(a)-10(c). Graphically, the operation (3.5) corresponds to flipping the diagrams of Fig. 
10 around the horizontal axis and changing the designations + ^ — on the fermion lines. 

The next step is to choose combinations of the three frequency integration variables that 
make the integrals as simple as possible. For each of the terms in Fig. 10, one can find 
three independent loops that have all arrows but one going around the loop in the same 
direction. Choose the frequency of the single oppositely oriented line of each such loop to 
be an integration variable. In Fig. 10, these three frequencies are shown explicitly for each 
term in the first row, where /q = —{po+qo+ko), uq = —{qo+ko), and vq = po+ko- By closing 
the frequency integration contours in the appropriate half plane, one can pick up a single 
pole for each corresponding to the labeled lines in the figure. For instance, from Fig. 10(a), 



^ Another such symmetry is ^ —{P+Q+K) with P and Q unchanged. 
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FIG. 10: The six nonzero terms of Fig. 6(a) arising from the decomposition (3.4) represented by 
Fig. 8. The conventions for defining the direction of frequency flow are not given by the arrows 
here for the fermion hnes but by the original diagram of Fig. 6(a) for pQ, qq, and ko and by the 
definition Iq = — (po+'Zo+^o) for Iq. 



one obtains the term 



dpo dqo dlo 



[po + Ep-ie) ^ 



pgk2Ep2Eq2Ek2Ep+g+kJ 2m 2m 2m 

X (-fco - Ek + i£y^{-lo -po-ko + Eg- iey^{-lo - Ep+q+u + 

X (Po + kQ- |ep+fc + i^Y'^iPo + lo- |eq+fc + i^)'^, (3.6) 

which integrates to 

Jpqk '^^p '^^q '^^k '^^p+q+k 

X [Ep + i^fc + \ep+k) ^{Ep + Ep^q^k + l^q+fc) ^- 



(3.7) 



Doing all the terms of Fig. 10 similarly, we obtain the following result for the contribution 
of the cross diagram to the effective potential: 



9^0 



{Ep 



pqk 1 2Eq 2Ek 2Ei Spi 



[P ^ Q) 



{Eq + tq) (Eq 



pk^pqkl 



^ pq.k 



1 



s., 



pk 



qk 



(3.8) 



where we introduce short-hand notation which will also be convenient for other diagrams: 

l = -{p + q + k), (3.9a) 
Spk = Ep + Ek + ^Cp+fc, (3.9b) 



■^pq,k 



E. 



Eq ~\~ 2^P+k ~l~ 2^9+'"' 



pqkl 



Ep + En + Ek + El. 



(3.9c) 
(3.9d) 
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The {p ^ q) term at the end of Eq. (3.8) represents the second row of Fig. 10 [using 
the {P,Q,K) —{Q,P,K) symmetry discussed earher combined with parity {p,q,k) — > 
— (p, q, k)]. One may drop the (p q) in favor of multiplying the rest of the expression by 
a factor of 2. 

Some readers, used to perturbation theory about free Fermi gases, may wonder at the 
absence of Fermi-sea step functions such as 9{ix — ep) in the expression (3.8). This is because 
of the effects of the condensate which is large compared to the chemical potential fi ~ e0o- 
The system is not describable as a small perturbation to a free Fermi sea. In Appendix B, 
we give a brief, illustrative example of what happens to a standard — e^) factor if one 
adds a condensate 0o and increases it to 0o ^ A*- 



2. The momentum integrals 

The remaining integrals (3.8) that we need to do can be made dimensionless by rescaling 
momenta as 

p^(2m(/)o)'/V, (3.10) 

which has the effect of replacing g^^^Q by g'^(j)Q^{2m(j)o)^^^'^ outside the integral and replacing 
Ep and Ep by the dimensionless versions 

ep^p', Ep^ip' + iy/' (3.11) 

everywhere inside the integral. Using the formula (2.5) for g'^, this can be written as 

Vr = <Po(^y\'JCr\e), (3.12) 

where 

X:(®)(e) = 4(47r)3'^/2 f '^{Ep-ep)^ (fi + ^g) iEq_ - eg) / I , M 

Jpqk'^Ep2Eq2Ek2EiSpi iSpkSpqki %pq,k \^pk ^qk J 

(3.13) 

is a dimensionless function of e. For a NNLO evaluation of the potential, we may evaluate 
JCf^ at e = 0, so that 

^r = 0o(^)'''eM^r + O(^)]> (3-14) 

where the numerical constant Kg®^ is the result of the integral (3.13) evaluated in exactly 
four dimensions. 

For numerical evaluation, we use rotational invariance to rewrite the momentum integrals 
as a six- dimensional integral over the magnitudes of and angles between the momenta: 

(2^ (2^ (2^ " (2^ Jo ^^'^'^'^^'^^^^ / sm^epsm^egsm^gdepde^d^g, (3.15) 
corresponding to a choice of Cartesian coordinates aligned so that 

fc = A;(1,0,0,0), (3.16a) 
p — p (cos ^p, sin ^p, 0, 0) , (3. 16b) 

q — g(cos^q,sin^qCOS^q, sin^qSin^q, 0). (3.16c) 
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We then perform the integral (3.13) numericaUy using adaptive Monte Carlo integration. 
The result is 

irf^ ~ 0.15101, (3.17) 

where ~ in this paper will mean that the result has an estimated error of at most ±1 in 
the last digit. Equations (3.14) and (3.17) represent our final result for the cross diagram of 
Fig. 6(a). 



C. Basic two- loop diagram through O(e^) 

The basic two-loop diagram of Fig. 5(a) is given by 

-iT/W = -{igf [ iG^iiP) iG22i-Q) iD{P + Q). (3.18) 

JPQ 



IPQ 

Doing the frequency integrals, 

Ko) _ f i^P - ^p){Eq - eg) 

Ipq '^Ep2EqSpq 

Rescaling momenta as in Eq. (3.10), 

d/2 



J„„ Zrj„ Zrjn On 



y^'-^o{^j elC^ie) (3.20) 
with 

JC^ie) = -2iAnr [ ■ (3.21) 

J pq ZHip ZHiq Opq 



If one sets d—A, one has /C2(0) = — C2 where 

C2 ^ 2{4n)' y ^ ^ , (3.22) 



d^p d^q {Ep - ep){Eq - e q) 
(2^ (2^ 2Ep2EqSpq 



which gives 

C2 ~ 0.14424 (3.23) 

upon numerical integration, all just as in Ref. [6]. For a NNLO calculation of the potential, 
however, we need the next term in the expansion of lC2{e) in e. We can obtain this by 
rewriting the d dimensional integration as 

7^ 7^ = TT- / P^-^q^'^dp dq / sin'^-^ 0pq d0pq, (3.24) 

(27r)'^ (27r)^ (4vr)'^+ir(f ) r(£^) Jo Jo ^ ^ 

where 9pq is the angle between p and q, and then expanding in e. The result is 

/C2(e) = - [1 + (I - 7 - In 2)e] + eCf + 0{e^) (3.25) 

where 

Ci°'^ = 2{4nr / 7^ 7^ iEp-ep){Eq-eq) ^.^^^ . ^ ^ ^ 

' ^ W (27r)4 (27r)4 2Ep2EqSpq ^ ' 

with numerical value 

C?°^^ ~ 0.14238 . (3.26b) 
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D. Scalar loop with two self-energy insertions 



The next diagrams we consider are those of Figs. 6(f)-6(h), which correspond to a scalar 
loop with two self-energy insertions U and the corresponding counter-term diagrams. In- 
dividually, the diagrams of Figs. 6(f) and 6(g) are ultraviolet (UV) divergent; it is only in 
their combination that the divergences are eliminated. So we must be careful not to set d=4 
in integrations until we have organized the terms into absolutely convergent integrals. We 
will keep d general in what follows until near the end. 

Together, Figs. 6(f)-6(h) give 



(nn) 



V 



[tD{V)]\-t[U{V)-flo{V)]}', 



(3.27) 



where V — {vo, v) is the frequency and momentum of the scalar line, U is the one- loop scalar 
self-energy given by 



-m{V) = -{zgf / zGn{P)iG22{P-V), 



(3.28) 



and the corresponding counterterm Ho is given by Eq. (2.7). 

The frequency integral in Eq. (3.28) is straightforward, especially using Eq. (3.4), yielding 



m) [ 
jp 



2£/p '2Ep—^ 



i^Ep -\- €p^ (^Ep—D "f- 6p_.i;) (-^p ^p){^Ep—V ^p—v) 



(3.29) 



vq- Bp - Ep_v + is vq + Ep + Ep_v - i£ 

The momentum integral would be UV divergent if we set d=4. We will isolate this divergence 
by isolating the large p behavior of the integrand by writing 

n(y) = ndiv(v) + n,eg(v), (3.30) 

where the divergent piece of the momentum integral is 

1 



ndiv(^) = 9 



and the remainder is 



pVo- Sp- Sp-v + ^£ 

i^Ep -\- €p) {Ep—i, -\- 6p— i;) (-^p ^p) {.-Ep—v ^p—v) 



(3.31) 



Vq - Ep- Ep^v is 

1 



vo + Ep + Ep^v — 



Vo — ^p — ^p-v + i£ 



(3.32) 



The strategy here is to have chosen the form of Ildiv to be simple enough that we can 
manage to evaluate the integral in general dimension d. The result of evaluating Eq. (3.31) 
is^ 



n, 



div 



(y) = -/r(i-f) (- 



-Vq + lev-i 



200 



-e/2 



(3.33) 



The dependence of (3.33) only on the combination —vq + ^e^ is a result of Galilean invariance. Galilean 
invariance is broken by the condensate ^o, but does not appear in our definition of Ildiv 
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where Eq. (2.5) has been used for g'^. Now use our spht (3.30) of H to rewrite the contribution 
to the potential given by Eq. (3.27) as 

= -i j|^[i^(i^)]'|[nreg(i^)]V2n.eg(i^)[ndiv(i^) -no(^)] 

+ [ndiv(V^)-flo(V^)]'|. (3.34) 

In dimensional regularization, the terms which do not involve Hreg integrate to zero because 
of their scaling properties. For example, changing integration variables {vq, v) — > (A^i'o, \v) 
for an arbitrary constant A, 

= / '^'^[^-'D(V)f[nm' (3.35) 

So this integral must equal itself times A'^*^'^^^^ and so must vanish. For similar reasons, we 
may dispense which each term generated by expanding the square in 

For the first term in Eq. (3.34), we continue by performing the Vq integration, using Eq. 
(3.32) for Ilrcg. One may avoid the bother of dealing with the double pole from D[K) by 
closing in the upper-half plane. The result is 



^y|^[D(v)]^[n,eg(v^)]' 

2 Jpqk 1 2£'p 2£'q 2Ek 2Ei S^^, 



Sql ^pqkl 

2Eq2Ei 



{Ep + €q + Ek + €l) 



+{pk^ql)}, (3.37) 



where we have used the notation of Eqs. (3.9) and the momentum naming conventions shown 
in Fig. 6(f) {so v = p + k). This integral is absolutely convergent in (i=4, and so we may 
evaluate it numerically just as we did for the cross diagram in Sec. IIIB2, giving 

-1 l[D(V)nn„,(V)r = ^„ [A-r" +0(.)] (3.38) 

with 

^^nn,i) _ 0.006753 . (3.39) 

Finally, we need the second term from the right-hand side of Eq. (3.34). Taking Ilreg and 
Ildiv from Eqs. (3.32) and (3.33), and closing the vo integration in the upper half plane to 
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avoid the branch cut in Hdiv, one finds 



[D{v)Yu,,^{v) [ndiv(^) - rio(^)] 

iE,-e,)iE,-e,) [n,,,^^) _ n,(y)] 



pfe '^Ep 2Ek 5^^, 
2 / {Ep — ep){Ek — ek) 



pk '^Ep 2Ek Spk 



ir(i-f) 



-e/2 

- 1 



(3.40) 



If we expand about four dimensions, we obtain absolutely convergent integrals at every order 
in e. We are therefore free to expand the integrand in e to obtain 



[l + 0(e)] 



eg' 



(Ep — €p)(Ek — efc) 



2 Jpk '^Ep 2Ek Spk 
Rescaling momenta in the usual way then gives 



1 — 7 — In 



S. 



pk 



200 



(3.41) 



V 



[D{V)fu,,g{V)[Uai.{V)-flo{V)] 



— (^0 



27r 



d/2 



{ [i(l - 7 + In 2)^2 - Ct''%' + 0(e3)} (3.42) 



with C2 given by Eqs. (3.22) and (3.23) and 



a 



(log,5) 



2(47r)' 



(iV d'^q {Ep - ep) {Eg - e,) ^ 



(27r)4 {2ny 



2Ep 2Eq Spq 



\n{Spq) 



(3.43a) 



givmg 



(log,5) 



0.19408. 



.2 (3.43b) 

Our final result for the contribution of Figs. 6(f)-6(h) to the effective potential is the sum 
of Eqs. (3.38) and (3.42): 

= 00 (^) { [K^''^ + 1(1 - 7 + In 2)C2 - C?-'^)] + 0(6^) }. (3.44) 



E. Remaining diagrams 

The remaining diagrams are relatively easy, and we summarize results in Appendix A. 
Here, we will just make a few comments on method. 

A simple way to handle diagrams with self-energy insertions is to consider the diagram 
without any such insertions, and then replace ep by = ep — in both ep and Ep = 
{€p + 0o)^^^ fermion energies, and replace |e^, by ^e^ — 2/x for boson energies. For 
instance, to simultaneously evaluate all of the two-loop diagrams of Fig. 5(b)-5(d) with 
chemical potential insertions, replace Eq. (3.19) for the basic two-loop result V-j^^ by 

_ 2 / (^-^^)(^^-^^) .3 45^ 

^ Jpq 2Ej; 2E!^ {Ef^ + E^ + iep+q - 2/x) ' ^ ' ^ 
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where E!^ = [(e^ — /x)^ + ^ol^^^. Then Taylor expand the integrand to the desired order in 
H, which in this case is first order. The same method can be used on the one-loop integrals 
of Fig. 4 starting from the basic one-loop integral of Fig. 4(a), ^ 

=i J lndet[-iG'-^(P)] = -J Ep. (3.46) 



IV. RESULT FOR ^ 



Combining the results given previously and in Appendix A, the full effective potential at 
NNLO in e is 



y(0o) = 



m4>o 
27r 



d/2 



^(1 + a,e + a,e') - ^(1 + b,e + b,e') - ^ 
6 e Z0O 



+ 0{e% (4.1) 



where 0o is treated as 0(1) and as 0(e). The various numerical coefficients are 
Qi = 1 (I - 7 - In 2) - 302 ~ 0.09877, 



02 



8 V3 



- 7 - ln2)' + if + 3 ( - 1 + ^7 + I ln2)02 + Of - O 



-((log,<S) 



(4.2) 

(4.3) 
(4.4) 
(4.5) 

(4.6) 

The ratio ^ can be computed from F(0o,A*) by the procedure used in Ref. [6]. First, we 
determine the expectation 0o which minimizes the potential: 



~ -0.15840, 

61 = |(|-7 + ln2) ^ 0.30797, 

62 = 1(1 - 7 + In 2)' + ^ - i^^'^) ~ 0.33703 , 



where 



-0.18348. 



= ^[l + 0.12586e + 0.56845 e2 + 0(e2)]. 



(4.7) 



Then we determine the fermion number density n from the pressure P — — F(0o); giving 

(4.8) 



n — 



dP dV 1 /m(^o^ "^'^ 



djj, dji e \ 271 



1 + 616 + 626' + ^ + 0(e=') 

00 



Next, we take the formula for the Fermi energy €^{71} of a d-dimensional ideal Fermi gas 
with the same density, 

(4.9) 



271 



^F = -[ir(| + i)n] 



2/d 



m 



There are similar issues as footnote 6 concerning the frequency integral and the behavior of its integrand 

at infinity. In general rcgularization schemes, this could be avoided by defining the integrals here with cf)Q 
(and /i) independent subtractions, which will not affect the determination of ^ from the effective potential: 
i /p{lndet[— iG^^(P; 0o)] — Indet [— iG'~^(P; 0)] } = — J^{Ep — ep). But, in dimensional rcgularization, 
the subtracted term vanishes anyway by scaling arguments similar to those reviewed in Sec. HID. 
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For the density (4.8), this then gives 



ep 00 



ir(i + i) 



l + 6ie + 62e' + ^ + C>(e3) 

00 



-2/d 



(4.10) 



Substituting in the expectation (4.7) for (po, and rewriting e^/'^ = e^/^'^e^/^, produces the final 
result (1.7) through NNLO. The additional logarithmic term shown in Eq. (1.7) at the order 
beyond NNLO is the subject of the next section. 



V. THE € EXPANSION BEYOND NEXT-TO-NEXT-TO-LEADING ORDER 
A. General 

In this section, we discuss a difficulty that arises in applying the diagrammatic expansion 
of Sec. II if one were to proceed to yet one higher order in e, attempting to evaluate the 
effective potential through O(e^) [and so ^ through 0(e^/^)]. The problem is an infrared 
problem arising from the fact that the scalar excitations, unlike the fermionic excitations, 
are not gapped. 

For the sake of specificity, consider the O(e^) contribution to the effective potential made 
by the diagram of Fig. 11(a). Together with the corresponding counterterm diagram, this 
gives a contribution to the effective potential 

^^(example) ^ _1 /" [D{V)YD{-V) [U{V) - Uo{V)] [ff(l^)]', (5.1) 

^ Jv 

where the one-loop self-energy 11 is given by Fig. 1 and 11 by Fig. 12. Let us now explore 
the contribution to this integral from small scalar frequency and momentum: |fo| ^ (f)o and 
<^ 00- At small V, we approximate V as zero inside the self-energies, and this region of 
integration contributes 

vt' ~ -n(o) [n(o)]'- / [d{v)Yd{-v), (5.2) 

^ Jv 

where the integration is restricted to small V. [Note that no(0) = 0.] The frequency integral 
is simple, giving 

vt' - n(o) [n(o)] ' /" 4 ■ (5.3) 

J V 

The momentum integral is IR divergent in d < 4. So our evaluation of diagrams has broken 
down in the infrared if 11(0) and 11(0) are nonzero. 



If one blindly tried to regulate the IR divergences using dinicinsional regularization (which we have previ- 
ously used only to regulate the UV), then the IR momentum integral in Eq. (5.3) would generate a factor 
of 1/e, which would in any case destroy the e counting of Sec. II. The integral does not generate zero 
(by scaling arguments like those of Sec. HID) because the 1/eJ integrand is an approximation valid only 
for \v\ <C 00- The behavior of the full result changes for > (bo, and so (po provides a scale. (We say 
"blindly tried to regulate" because dimensional regularization throws away nonlogarithmic divergences, 
and it does not distinguish between IR and UV logarithmic divergences. This means it can be a dangerous 
procedure unless you already know that all divergences will cancel in the final result.) 



18 



*' .-X- 

FIG. 11: Two examples of diagrams that produce logarithmic infrared divergences. 




FIG. 12: The self-energy 11 that mixes (p with (/?*. 



Another example of a diagram producing similar problems is shown in Fig. 11(b). Adding 
one or more additional self-energy or chemical potential insertions to the scalar loop in either 
diagram would generate even more severe (power-law) infrared divergences. 

We shall see in a moment that 11(0) and 11(0) are both 0(/x). The problem, then, is that 
when V is small enough that e„ is of order <^ 0o, it is no longer a good approximation 
to treat the scalar self-energies, or the chemical potential 2//, as perturbations to — 
—Vq + One must therefore rcsum 11(0), 11(0), and 2/x into scalar propagators in order 
to recover a well-behaved perturbation theory. Equivalcntly, one must use the appropriate 
low-energy scalar effective theory when evaluating the contribution of low-energy scalars 
{\vq\ <S 0o) to the effective potential. At low energy, the effective scalar propagator would 
be generated by the effective interactions 

£eff - [idt + ^) V + 2^xip*Lp - n(o) LP* LP - n(o) \{lplp + (^v*) 

with corresponding propagator 

^^p,(p^-\ep + 2jx-Il{Qi)+ie -n(0) 

^ ' V -n(o) -po - \ep + 2/1- n(o) + iej ■ ^ ■ ^ 

To better understand the structure of this propagator, we now evaluate 11(0) and 11(0). 
Prom Eq. (3.29), 

n(0) is given by 
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In dimensional regularization, the momentum integrals give 



and 



n(o) = i^n(o) + 0(e)], 



(5.8) 



(5.9) 



{d + 2) 

where we have used Eqs. (2.5) and (4.7) for g"^ and 0o- To leading order in e for scalars with 
energy of order /i, the effective low-energy interactions (5.4) are then 



with 



V{P) 



11 + ie 



2 P 



H + ie 



-1 



(5.10) 



(5.11) 



The imaginary part of (p remains gapless, which must occur in any consistent approximation 
since it corresponds to a Goldstone boson. 

With this formula, we can now compute the leading contribution to the effective potential 
from low-energy scalars and see that all is well. Analogous to the fermionic result (3.46), it 
is 

y(soft^)^_^ I lndet[-iV-\P)] ^ J ln[pl - ^e^C^e^ + 2t,) + is] 

-1 1/2 

|ep(|ep + 2//)J . (5.12) 



This integral has no infrared divergences. If we naively expanded the integrand in powers of 
fi, wc would obtain infrared divergences starting at third order ~ /i^ Ip^p^- This is the same 
order as the result (5.3) for Fig. 11 that started this discussion. If we instead integrate Eq. 
(5.12) up to a UV momentum cutoff (4mA) for the effective theory, we find an expansion 
in n of the form 



yisoft (fi) 



m 



#A3 + + #/^'A + #/^' In (^) + #/^' 



+ 



(5.13) 



6o, so one should very roughly think of 



The soft effective scalar theory breaks down at ep > 
the energy scale A as of order 0o- 

The moral of this story is that one will have to make a proper treatment of low-energy 
scalar fields in order to go to higher orders in e than the NNLO calculation performed 
in the bulk of this paper. In particular, the approach of Sec. II would lead to logarith- 
mic divergences in the effective potential at O(e^) and worse divergences at higher or- 
der, but a proper resummation of soft scalar physics will resolve these divergences into 
0(/x^ ln(/x/0o)) = O(e^lne). This corresponds to a correction to ^ of order e^/^'^e^/^ In e. 
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B. The coefficient of the logarithm 



Though we are not prepared to make a full calculation of the order beyond NNLO in 
the e expansion, it is simple to extract the coefficient of the logarithm at that order from 
the preceding discussion. Naively expand the integrand in the formula (5.12) for the soft 
contribution v') third order in fj,. The third-order term is 



(5.14) 



In four dimensions, this is the logarithmically divergent integral 



87r2 



dp 



2 S 

m jj, 
47r2 



de. 



(5.15) 



But, from the discussion surrounding Eq. (5.13), we know that this logarithm is cut off in the 
infrared by the energy scale due to the necessity of resummation, and in the ultraviolet by 

the energy scale 0O) the energy scale where the self-energies are no longer well approximated 
by their zero-momentum values. Thus, even though we cannot easily compute the constant 
under the logarithm, we can write^^ 



9 

m n 
47r2 



Inl^KOdl 
/i 



(5.16) 



We can now use this to obtain the explicit logarithm shown in our final result (1.7) for ^ 
beyond NNLO by including Eq. (5.16) in the analysis of Sec. IV. Because dV^^ V^0o does 
not have a logarithm, there is no logarithm in the NNNLO result for the location 0o of the 
minimum of the effective potential. The logarithm appears in ^ only through its effect on 
the density n — —dV/d/i. Equation (4.8) for n is modified to 



n — 



1 [ m(j)Q 
2^ 



d/2 r 



1 + 6ie + 626^ + 



ln^ + 0(e3) 



(5.17) 



givmg 



ep 00 



r(| + i) 



1 ^ ^ 2 , e/i 3e/x2 00 , 3^ 
1 + 016 + 026 +- r^ln \-0(e^ 



00 00 



~2/d 



(5.18) 



Substituting the expectation (4.7) for 0o then produces the NNNLO logarithm shown in the 
final result (1.7) for ^. 



VI. EXTRAPOLATION OF ^ AT d=3 

Because of the large relative size of the 0(e^/^) term in our result (1.7) for ^, the next- 
to-next-to-leading order (NNLO) result will likely only be useful in conjunction with more 
sophisticated analysis of dimension dependence than the naive prescription of setting e—1. 



The coefficient of this logarithm has also been computed by Y. Nishida [17]. 
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In this respect, the situation appears somewhat analogous to the e expansion for the critical 
exponent uo in the Ising model, u is the exponent characterizing corrections to scaling, and 
its e expansion is [12] 

c^ising = e - 0.62963 + 1.61822 + O(e^). (6.1) 

Despite the large size of the NNLO term, historical analysis using the terms shown 
above gave u ~ 0.79 with a simple Borel-Pade approximation [9]. The latter is within a few 

percent of the correct result. 

Unfortunately, the ratio ^ docs not have as simple an analytic structure as do critical 
exponents. Our result for ^ may be written as 

e=ie^+2/'^6(e) (6.2) 

with 

6(e) = 1 - 0.04916e - 0.95961 - | Ine + 0{e^''^) (6.3) 

Critical exponents have a simple asymptotic expansion in powers of e. But, even if we factor 
out the overall e^+^Z"^ and focus only on 6(e) above, our expansion contains powers of Ine. 
The Borel transform of a series 

f{e) = yjne- (6.4) 



is the faster-converging series 



F(t) = J]^r. (6.5) 

n 

The original /(e) may be recovered from its Borel transform by 

poo 

/(e) = dt e-*F{et) (6.6) 
Jo 

if F{t) does not have any singularities on the positive real axis which make the integral 
ill defined. The standard approach for critical exponents is to fit some type of Pade-like 
approximation to the Borel transform. A simple [M/N] Pade approximation would be 



where we assume /(e) is normalized so that /(O) = 1. More sophisticated versions have been 
used for accurate estimates to critical exponents from high-order e expansions. However, 
such estimates must necessarily break down for the Borel transform B{t) of (6.3) because 
of the Ine terms. The appearance of Ine in the small e expansion of 6(e) gives rise to Int 



See, for instance, the results for 9 = ivu and v from three-dimensfonal (3D) series techniques, the e 
expansion, numerical Monte Carlo simulations, and experiments all reviewed in Ref. [12] for the 0{N) 
model. (The Ising universality class corresponds to = 1 in these tables.) The current results presented 
there for co from 3D series techniques and the e expansion are 0.799(11) and 0.814(18). 
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FIG. 13: Extrapolations of ^ vs dimension d. The thin sohd hnes show the result of Pade-Borel 
extrapolations of type [4/0], [2/2], [3/1], and [0/4] from bottom to top. For comparison, the thick 
solid lines show simple truncations of the expansions about four (right) and two (left) dimensions, 
as discussed in the text. The thick dashed line is the truncated d = 4 — e expansion at NLO rather 
than NNLO. 

terms in the small t expansion of its Borel transform, and these are not accommodated by 
Pade approximants such as (6.7).^^ 

Clearly, what is needed to make full use of e expansion results is a full understanding of 
the analytic structure of ^ (e) in e, in order to inform the strategy for how best to extrapolate 
to (i=3. Nonetheless, it is interesting to see what happens if one naively extrapolates ^(e) 
using simple Pade estimates (6.7) to the Borel transform B{i). This was carried out by 
Nishida and Son in Ref. [13], constraining the Pade approximants by (i) the next-to- leading- 
order (NLO) result near the expansion of ^ in 4— and (ii) a next-to- leading- order result 
for the expansion of ^ in c? — 2. Since the leading-order 4— (i result has already been used 
to normalize 6(0) = 1 above, this represents three constraints on -B(t), and so the possible 
Pade approximants are those with M -\- N — 3. Extrapolating to three dimensions, they 
then found i = 0.391, 0.364, and 0.378 with [3/0], [1/2], and [0/3] Pade approximants. They 
did not find any solution satisfying the constraints for a [2/1] approximant. These values 
for i span 0.378 ± 0.014. 

If we naively follow the same procedure but add the information from our NNLO coeffi- 
cient, we find I = 0.300, 0.367, 0.359, and 0.376 in three dimensions from [4/0], [3/1], [2/2], 
and [0/4] approximants. We did not find a solution for [1/3]. The [4/0] value is an out-lier, 
which makes a certain amount of sense. [4/0] corresponds to a simple polynomial form for 



Specifically, the Borel transform of 6(e) = 1 + ae + /Se^ + 76^ In e + (5e^ \s B{t) = l + at+ ^(3t^ + ^jt^\nt + 
^[i5 — 7 V'(4)]<^, where V(-z) is the digamma function. This can be demonstrated by checking Eq. (6.6). 
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FIG. 14: A selection of estimates of ^ in tliree spatial dimensions. Experimental values (circles) 
include 0.74(7) [18], 0.34(15) [19], 0.321°;}^ [20], and the more recent 0.51(4) [8] and 0.46(5) [21]. 
Simulations include fixed-node Green's function and diffusion Monte Garlo upper bounds of 0.42(1) 
[4] (square). Other simulation methods have estimated 0.22(3) [22], 0.07 < ^ < 0.42 [23], and ^ 0.44 
[24] (the last two are based on finite-temperature calculations). The NLO 4 — e value is that of 
Eq. (1.6) [6]. The rest are the Pade-Borel estimates discussed in this section, with the purely NLO 
result from Ref. [13]. 

B{t). This does not allow for any singularities in the Borel plane and so will not produce 
an asymptotic series in e, whereas the large NNLO coefficient suggests that the asymptotic 
nature of the expansion should not be ignored at this order. If we focus only on the other 
values, they span ^ = 0.367 ± 0.009. This is consistent with the NLO results, but not an 
obvious improvement. The fits of ^ as a function of d are shown in Fig. 13 and are quite 
similar to the earfier fits of Ref. [13] discussed above which did not use the NNLO result for 
d = A-e. 

For comparison, the same figure also shows the result of avoiding any fancy extrapolation 
but simply naively using the truncated NNLO result C = |e^+^/'^[l - 0.04916 e - 0.95961 e^] 
for d = 4 — e or the corresponding NLO result ^ = 1 — e of Ref. [13] for d = 2 + e. It is 
amusing that both of these naive extrapolations happen to give consistent values of ^ ~ 
at d = 3. We imagine that this is a coincidence. 

Figure 14 summarizes a selection of estimates of ^ from experiment and numerical sim- 
ulations, and compares them to results obtained so far from the e expansion. The current 
Pade-Borel extrapolated e expansion results, however, need to be taken with a significant 
grain of salt because their assumptions of analytic structure are not consistent with what 
we have learned in Sec. V about higher-order corrections. Further study will be required to 
develop more consistent extrapolations. 
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APPENDIX A: SUMMARY OF RESULTS BY DIAGRAM TYPE 



Fig. 4(a): v!"> = - / = ^„ (!^)' 



,,0, , , . i'm,p„Y''' n-i-i) 



+ [1(|-T-In2f + l£VO(£»)} (Al) 
+ [i(i-T + ln2)'' + i]£2 + 0(a} (A2) 



Fig. 5(a): See Eqs. (3.19) through (3.23). 

Figs.5(b-d): Vi^^ ^-^^g^ll (^ + 



pq 2-E'p ^pq \ Eq ^pq 



W^) iff +0(^1 (A4) 



d/2 

. 27r '^ 

K^^''^ ~ -0.25835 (A5) 
Fig. 6(a): See Eqs. (3.8), (3.14), and (3.17). 
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Figs. 6(b-d): V, 



Fig. 6(e): 
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(A8) 
(A9) 

Figs. 6(f-h): See Eqs. (3.44), (3.39) and (3.43). 

We have used the shorthand notation of Eq. (3.9). In the momentum integrals shown in 
Eqs. (A4) and (A6), we have used symmetry to simphfy the expressions. For expressions 
that correspond to simply doing the frequency integrals of these diagrams without using 
such symmetry, simply symmetrize the integrands under p ^ q. 



APPENDIX B: WHAT HAPPENED TO 6>(At - ep)? 

In this appendix, we give a quick, illustrative example of the effect of 0o on what oth- 
erwise would be a Fermi-sea step-function 6{fi — ep). The example will be the mean-field 
theory result for the fermion number density n in the presence of both a non-negligible 
chemical potential /j, and a condensate 0o- In terms of diagrams, this corresponds to the 
one-loop fermion diagram of Fig. 4(b) if we take the cross to represent the fermionic num- 
ber operator and (unlike the rest of this paper) include the chemical potential fi in the 
fermion propagator rather than treating it as a perturbation. It will be clearer if we give 
a regularization-independent version of the result rather than continuing to work in dimen- 
sional regularization. In terms of diagram evaluation, this can be achieved by subtracting 
the vacuum contribution {/i — (f) — 0) to n before doing any integrals. The result, after 
frequency integration, is then the standard mean-field formula for the number equation^^ 

In the limit 0o 0, the above formula gives the usual result for a free fermion gas: 



n 



2 f 9{ii-ep). (B2) 
Jp 



^'^ For comparison, the result to leading order in /i in dimensional regularization could be obtained by 
applying n — —dV/dji to the result Vi = — J^E^ discussed in section HIE, giving n = — J^et^/ in 
dimensional regularization. This appears to differ from Eq. (Bl) by 1, but 1 vanishes in dimensional 
regularization. 
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It would be problematical to treat /i perturbatively in such an expression. However, in 
the opposite limit ji <^ relevant to this paper, there is no obstruction to treating 
perturbatively. To leading order in /x. 




We thus find fractions involving energies rather than d functions, similar to results we have 
derived for other diagrams in this paper. 
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